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1 .  INTRODUCTION 


Most  existing  schemes for  computing  transonic  potential  flowfields 
about  wings  or  wing-body  combinations  are  limited  to  first-order  accuracy  in 
supersonic  regions  because  of  the  use  of  first-order  upwind  differencing  or 
the  addition  of  first-order  artificial  viscosities  or  densities.  Although 
these  schemes  have  been  used  extensively  for  preliminary  aerodynamic  design, 
most  of  their  successful  applications  are  limited  to  simple  wing-body  geom¬ 
etries.  Extension  to  more  complex  geometries  is  prohibited  by  the  need  for  a 
large  number  of  grid  points  for  adequate  prediction  of  shock  location  and 
strength,  especially  in  the  case  of  double  shocks  which  are  common  features  of 
flow  about  highly  swept  and  aft-cambered  wings.  Development  of  higher-order 
schemes  is  therefore  necessary  to  improve  the  resolution  of  solutions,  with  a 
reasonable  number  of  grid  points,  for  realistic  wing-body  geometries. 

Several  second-order  schemes®"^  have  been  introduced  for  airfoil  and 
cascade  flowfield  calculations.  Jameson*-®  and  Chen^  demonstrated  that  second- 
order  fully  conservative  and  quasi-conservative  schemes,  respectively,  are 
capable  of  predicting  double  shocks  on  an  airfoil  surface,  which  cannot  be 
accurately  resolved  using  a  first-order  scheme  without  a  large  number  of  grid 
points.  Chen^  also  demonstrated  that  his  second-order  quasi-conservative 
scheme  provides  better  resolution  of  a  double  shock  than  the  second-order 
fully  conservative  scheme.  Ives  and  Liutermoza®  showed  that  their  second- 
order  nonconservative  scheme  provides  better  resolution  for  transonic  cascade 
flows  than  first-order  nonconservative  schemes.  A  discussion  of  first-  and 
second-order  nonconservative  schemes  has  also  been  given  in  Reference  7.  A 
study  of  artificial  viscosities  and  conservative  shock-point  operators  of 
different  orders  was  provided  by  Chen  and  Caughey,**  who  also  introduced  a 
third-order  quasi-conservative  scheme.  In  the  present  study,  second-  and 
third-order  artificial  viscosities  are  first  Introduced  for  transonic 
potential  flowfield  computations  about  wings  and  wing-body  combinations. 

Methods  for  differencing  the  small  disturbance  equation  at  shocks  were 
Investigated  by  Murman*^  and  Hafez.*®  Methods  for  treating  shocks  in  a  full 
potential  formulation  were  studied  by  Jameson*^  and  Chen  and  Caughey.**-  Fully 
conservative  schemes  for  treating  the  potential  equation  conserve  mass  flux 


isent Topically  across  shocks;  therefore,  the  predicted  shocks  are  always 
stronger  than  Rankine-Hugoniot  shocks. ^  In  the  present  study,  a  shock-point 
operator  is  derived  from  an  approximate  one-dimensional  flow  analysis.  Use  of 
this  operator  results  in  Mach  number  jumps  at  shocks  which  are  in  reasonable 
agreement  with  Rankine-Hugoniot  values.  Methods  for  differencing  at  shocks 
are  evaluated  using  a  model  problem  consisting  of  flow  through  a  converging- 
diverging  planar  channel,  with  a  shock  in  the  diverging  section.  In  flowfield 
calculations  about  wings,  partially  conservative  shock-point  operators  provide 
results  that  are  in  better  agreement  with  experiment  than  a  conservative 
scheme.  To  the  author's  knowledge,  this  result  is  1)  the  first  demonstration 
of  a  successful  third-order  scheme  applied  to  solution  of  the  full  potential 
equation,  2)  first  presentation  of  both  second-  and  third-order  solutions  for 
transonic  potential  flowfield  computations  about  wings  and  wing-body  configu¬ 
rations,  and  3)  first  demonstration  of  a  shock-point  operator  that  produces 
Mach  number  jumps  in  a  potential  flow  which  are  in  reasonable  agreement  with 
Rankine-Hugoniot  values . 
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2.  NOMENCLATURE 


speed  of  sound 
stagnation  speed  of  sound 
speed  of  sound  at  H  ■  1 

transformation  matrices  defined  in  Equations  (3),  (4),  and  (5) 
coefficients  defined  in  Equations  (9)-(17)  and  also  Equations 
(38)-(42) 

determinant  of  Jacobian  transformation  matrix  defined  in 
Equations  (33)  and  (47) 

reduced  velocity  potential  defined  in  Equation  (69) 
coefficients  defined  in  Equations  (21)-(29) 
artificial  viscosities  at  supersonic  points  defined  in 
Equations  (53)  and  (56) 

artificial  viscosities  at  shock  points  defined  in 

Equations  (55)  and  (57) 

mass  flow  rate  defined  in  Equation  (68) 

local  Mach  number 

normalized  Mach  number  »  u/a * 

partially  conservative  parameter  defined  in  Equations 
(55)  and  (57) 

second-order  transformation  derivatives  defined  in  Equations 
(30)-(32)  and  (48)-(49) 
coordinate  in  streamwise  direction 

velocity  component  in  x,  y,  and  z  directions  defined  in 
Equations  (34)-(36) 

velocity  components  defined  in  Equations  (18)-(20) 

coordinates  in  physical  space 

location  of  shock  in  the  channel  axis 

switch  function  defined  in  Equation  (54) 

local  flow  density 

flow  density  at  M  >  1 


3.  FULL  POTENTIAL  EQUATION 


Quasi-conservative  schemes  are  used  to  solve  finite-difference  approxi¬ 
mations  of  the  full  potential  equation.  Therefore,  it  is  convenient  to  first 
formulate  the  full  potential  equation  in  computational  coordinates.  By  apply¬ 
ing  the  chain  rule,  derivatives  of  the  potential  function  <t>  in  physical  co¬ 
ordinates  (x,  y,  z)  can  be  related  to  its  derivatives  in  an  arbitrary  curvi¬ 
linear  coordinate  system  (X,  Y,  Z)  as  follows: 


where 


-1 


yx2 

2xxyx 

zx2 

2xxzx 

ryxzx 

A 

2xYyy 

zy2 

2ZyZy 

2yYzY 

yyyy 

xxyY+xYyx 

zxzy 

xXzY+xYzX 

yxzY+yy: 

yz2 

2xzyz 

Zz2 

2xzzZ 

2yzzz 

yxyz 

xxyz+xzyx 

zXzZ 

xXzZ+xZzX 

yxzz+y2; 

yxyz 

xyyz+xZyY 

zYzZ 

xYzZ+xZzY 

yYzz+yz: 

xXX 

XYY 

xXY 

xzz 

Xxz 

XYZ 

C  * 

yxx 

yYY 

yxx 

yzz 

yxz 

yYZ 

zxx 

ZYY 

ZXY 

Zzz 

Zxz 

zYZ 

The  full  potential  equation  to  be  solved  is 


(a2  -  u2)^  +  (a2  -  v2)^  +  a2  -  w2)*zz  -  2uv<|>xy 


"  2vw*yz  -  2uw*xz  -  0  , 


where  u,  v,  w  are  the  x,  y,  z  components  of  the  flow  velocity,  respectively, 
and  a  Is  the  local  speed  of  sound  determined  from  the  energy  equation 

a2  -  a2  -  I—  —  (u2  +  v2  +  w2)  ,  C 

O  L 


where  y  is  the  ratio  of  specific  heats  for  the  assumed  calorically  perfect  gas 
and  aQ  is  the  stagnation  speed  of  sound. 

Substituting  Equations  (1)  and  (2)  into  Equation  (6),  and  after  perform¬ 
ing  matrix  inversion,  multiplication,  and  careful  algebraic  manipulation,  a 
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full  potential  equation  multiplied  by  the  determinant  of  the  Jacobian  trans' 
formation  matrix,  D,  in  arbitrary  curvilinear  coordinates  can  be  derived  as 


C1<),XX  +  C2^YY  +  C3^ZZ  +  C4^XY  +  C54>YZ  +  C6<I>XZ 
+  C7^X  +  C8^Y  +  C9^Z  =  ®  * 


where 

cx  -  [a2(h2  +  h2  +  h2)  -  U2]/D 

c2  -  [a2(h2  +  h2  +  h2)  -  V2]/D  ( 

c3  =  [a2(h2  +  h2  +  h2)  -  W2]/D  ( 

c4  -  [2a2(h1h/,  +  h2h5  +  h3h6)  -  2UVJ/D  ( 

c5  =•  [2a2(h4h?  +  h5hg  +  h6h9)  -  2VWJ/D  ( 

c6  «  [2a2(h1hy  +  h2hg  +  h3h9)  -  2UW]/D  ( 

c7  »  (h^  +  h2PY  +  h3pz^°  ( 

C8  “  ^h4pX  +  h5PY  +  h6PZ^°  ( 

C9  *  (h7PX  +  h8PY  +  h9PZ)/D’  < 

U,  V,  W  are  velocity  components  defined  as 

U  *>  hju  +  h2v  +  h^w  ( 

V  -  h^u  +  h5v  +  h6w  ( 

W  -  h?u  +  h8v  +  h9w,  ( 


coefficients  hj,  h2, 
defined  as 


h9  are  first-order  transformation  derivatives 


hi  »  yYzz  “  ^ZzY  (2^ 

^2  “  zYxZ  “  zZxY  (22) 

^3  =  XY^Z  “  XZ^Y  (23) 

h4  "  yzzx  ■  yxzz  (24> 

h5  -  zZxX  "  zXxZ  (25) 

h6  "  xz?x  "  xX?Z  (26) 

hy  ®  yxzY  “  yyzx  (27) 

hg  =  zXxY  —  zYxx  (28) 

hg  a  zxyY  —  xYyx  i  (29) 

and  coefficients  px,  pY,  and  pz  are  second-order  transformation  derivatives 
defined  as 

PX  C1XXX  +  C2XYY  +  C3XZZ  +  C4XXY  +  C5XYZ  +  C6XXZ 

PY  =  ClyXX  +  C2yYY  +  C3yZZ  +  C4yXY  +  C57YZ  +  C6yXZ  (31) 

PZ  “  C1ZXX  +  C2ZYY  +  C3ZZZ  +  C4ZXY  +  C5ZYZ  +  C6ZXZ  *  (32) 

The  determinant  of  the  Jacobian  transformation  matrix  is  defined  as 

D  -  hjX^  +  h4xY  +  h?xz  .  (33) 

The  velocity  components  u,  v,  w  are  defined  as 

U  "  <hl*X  +  h4^Y  +  h7^Z^°  O4) 

v  “  (^x  +  h5^Y  +  h8^Z^D  (35) 
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W  =  ^h3^X  +  h6^Y  +  h9^z^'°* 

Equation  (8)  can  be  reduced  to  a  two-dimensional  equation: 
C1(<,XX  +  C2^YY  +  C4^XY  +  C7^X  +  C8*Y  “  0  » 


where  Cj,  C2,  c^,  c-j,  and  Cg  reduce  to 
cx  -  [a2(x^  +  y2)  -  U2]/D2 

C2  "  ta2(xX  +  yX}  ”  y2 1/1)2 

c4  -  -  2[a2(x^xy  +  yxyY)  -  UV]/D2 

Cy  -  (xypY  -  yYPx^D 

c8  ■  <yXPX  -  Vy)/D 

and  U,  V,  u,  v,  J,  D,  px,  and  pY  are  redefined  as 
U  -  uyY  -  vxy 

V  -  uyx  -  ”x 

U  ■  (yY*X  -  yX*Y)/D 

V  -  (xXty  “  Xy<l>x)/D 

D  “  xXyY  “  XYyX 

PX  “  C1XXX  +  C2XYY  +  C3XXY 

PY  “  CiyXX  +  C2yYY  +  C3yXY  ’ 

Equation  (37)  is  consistent  with  the  two-dimensional  equation  derived  in 
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After  the  physical  coordinates  of  grid  points  have  been  prescribed,  the 

transformation  derivatives  Xj[,  *y»  XZ»  ^X*  *  *  xXX»  xXY*  XYY  *  *  *  can  *je 
computed  at  each  control  point  within  a  local  mesh  element.  A  second-order- 
accurate,  finite-difference  approximation  of  the  transformed  full  potential 
equation  thus  can  be  obtained  by  applying  a  second-order  element  (Figure  1). 
Within  the  element,  X,  Y,  and  Z  vary  from  -1  to  l  from  nodal  point  to  nodal 
point.  Therefore  the  mesh  element  is  uniformly  spaced  in  the  computational 
space.  Second-order  shape  functions  can  be  constructed  that  relate  the  func¬ 
tion  at  any  point  p  within  the  element  to  the  values  of  the  function  at  27 
nodal  points.  If  the  control  point  is  chosen  to  be  X  ■»  Y  *  Z  =  0,  then  the 
well-known  second-order,  centered,  finite-difference  formulations^*^  are  ob¬ 
tained  . 


Figure  1.  Transformation  of  a  second-order  element. 
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4.  LOCAL  SUPERSONIC  FLOW  AND  SHOCKS 


The  finite-difference  approximation  to  the  full  potential  equation  dis¬ 
cussed  thus  far  is  adequate  for  flows  that  are  entirely  subsonic.  To  treat 
transonic  flows,  proper  artificial  viscosities  or  densities  are  normally  added 
to  the  finite-difference  approximation  of  the  potential  equation  solved  at 
supersonic  points.  The  directional  bias  of  supersonic  flows  thus  can  be  re¬ 
flected  in  the  governing  equation. 

The  so-called  fully  conservative  and  quasi-conservative  schemes  conserve 
the  artificial  viscosities  or  densities  along  streamlines;  in  other  words,  the 
total  summation  of  artificial  viscosities  or  densities  added  to  the  potential 
equation  at  all  points  along  the  streamline  is  exactly  zero.  Naturally  the 
shocks  thus  predicted  are  consistent  with  lsentropic  mass-conserving  shocks, 
which  do  not  simultaneously  conserve  momentum  and  therefore  are  stronger  than 
the  Rankine-Hugoniot  shocks.  However,  in  the  nonconservative  schemes,  the 
total  summation  of  artificial  viscosities  or  densities  added  along  the  stream¬ 
line  is  not  zero.  In  the  limit  of  the  finite-difference  approximation,  this 
unbalanced  summation  term  appears  as  a  nonzero  source  term  on  the  right  side 
of  the  »vj*ential  equation  solved  at  certain  points  along  the  streamline  gen¬ 
erally  at  the  shock  point,  the  first  subsonic  point  downstream  of  the  shock. 
This  source  term  in  the  equation  represents  a  mass  source  in  the  flowfleld. 
Therefore,  the  solutions  thus  obtained  deviate  from  the  lsentropic  mass-con¬ 
serving  solutions. 

Since  the  potential  formulation  does  not  permit  simultaneous  conservation 
of  mass  and  momentum  flux  at  shocks,  errors  in  the  jumps  in  fluid  properties 
are  inevitable.  A  desirable  method  from  an  engineering  point  of  view  is  one 
in  which  errors  in  the  properties  of  primary  Interest  are  minimized.  As  a 
matter  of  fact,  it  has  been  consistently  shown  that  the  shocks  computed  by 
these  nonconservative  schemes  agree  better  with  experiment  than  those  computed 
by  the  conservative  schemes.^  The  nature  of  the  nonzero  source  term  was 
understood  to  be  related  to  the  addition  of  mass  flux;  however,  an  adequate 
mathematical  explanation  of  its  effect  has  not  been  given.  An  attempt  to 
explain  the  nonzero  source  term  will  be  given  in  Section  4.2,  in  the  form  of  a 
simple  one-dimensional  flow  analaysls,  following  the  introduction  of 
artificial  viscosities  in  Section  4.1. 
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4.1  Artificial  Viscosities  and  Partially  Conservative  Shock  Point  Operators 
The  second  derivative  of  the  potential  function  in  the  streamwise  direc¬ 
tion,  s,  is  given  as 

*88  "  ^2  (u2*xx  +  v\y  +  w\z  +  2uv*xy  +  2vw*yz  +  2uw<|,xz)  *  (50) 

Substituting  Equations  (1)  and  (2)  into  Equation  (50)  yields 

■  -T  (“2*xx  +  v2*n  +  “Sx  +  2UV*xy  +  2™*V2  +  20W»X2>  -  (51) 


where  U,  V,  W  are  given  in  Equations  ( 18)— ( 20) • 

The  directional  bias  of  supersonic  flows  can  be  properly  simulated  by 
performing  an  upwind  differencing  or  adding  artificial  viscosities  in  the 
approximate  streamwise  direction.  If  Y  ■  constant  lines  are  in  the  approxi¬ 
mate  s  direction,  the  principal  part  of  $  can  be  approximated  by 

s  s 


*ss  "  ~^Z  *XX  • 

A  second-order  artificial  viscosity  can  be  expressed  as 


r  ““^xx  1  /‘‘^♦xx  \  /““2»xx  \  /“"V\ 

-HH*  ■  i-n  * 


where 


p  -  max 


H  is  then  added  to  the  finite-difference  representation  of  Equation  (8)  at 
supersonic  points.  At  shock  points,  i.e.,  the  first  downstream  subsonic 
points  after  the  shocks,  the  following  first-order  artificial  viscosity  Hg  is 
added  with  pm  controlling  the  nonconservative  differencing: 


II 


(55) 


2 

If  pm  -  0,  the  quantity  pU  is  conserved  along  Y  *  constant  lines,  implying 
that  the  added  artificial  viscosities  are  conserved  along  approximate  stream¬ 
lines.  If  pm  >  0,  a  numerical  mass  flux  is  introduced  at  shocks,  modifying 
the  locations  and  strengths  of  the  shocks.  The  effect  of  pm  on  the  captured 
shocks  will  be  discussed  later.  Although  p  is  a  ramp  function,  both  H  and  Hs 
reduce  to  zero  as  the  mesh  size  goes  to  zero.  The  solution  is  second-order 
accurate  at  both  subsonic  and  supersonic  points,  and  first-order  accurate  at 
shock  points.  The  scheme  is  second-order  quasi-conservative.  In  the  so- 
called  quasi-conservative  schemes,  only  the  differencing  of  artificial  viscos¬ 
ities  is  in  divergence  form;  the  differencing  of  the  governing  potential 
equation  is  not.  A  second-order  fully  conservative  scheme  also  can  be  con¬ 
structed  by  incorporating  H  and  Hg  into  the  existing  finite-volume  algorithm. 

Third-order,  quasi-conservative  and  fully  conservative  schemes  can  be 
developed  by  adding  the  following  third-order  artificial  viscosity  at  super¬ 
sonic  points: 


and  adding  the  following  second-order  artificial  viscosity  at  shock  points: 


12 


If  pm  Is  set  to  zero,  the  quantity  is  conserved  along  Y  =*  constant 

lines.  If  pa  is  set  to  be  greater  than  zero,  a  numerical  mass  flux  is  added 
at  shocks,  as  in  the  second-order  schemes. 

4.2  A  Simple  One-Dimensional  Flow  Analysis  and  a  Shock-Point  Operator 

An  alternative  new  method  to  formulate  a  shock-point  operator  is  de¬ 
scribed  in  the  following  paragraphs.  This  method  is  based  on  an  approximte 
one-dimensional  analysis  in  which  information  from  the  Ranklne-Hugoniot  re¬ 
lations  is  Incorporated. 

The  flows  upstream  and  downstream  of  the  shock  can  be  considered  as  re¬ 
lating  to  two  branches  of  isentropic  flows.  Because  of  the  entropy  increase 
across  the  shock,  the  stagnation  density  decreases,  while  the  stagnation  speed 
of  sound  remains  unchanged  because  the  process  is  adiabatic.  The  continuity 
equation  can  therefore  be  written  as 

p*Mi(l  +X_=j.  MJ)"3  -  p^l  m2)“3  (58) 


or 


0, 


(59) 


where  x.  and  x.  are  the  axial  positions  just  upstream  and  downstream  of  the 
shock  and  p^,  p£>  and  M2  are  sonic  densities  and  Mach  numbers  of  the  flows 
just  upstream  and  downstream  of  the  shock.  If  the  Mach  number  is  assumed  to 
change  smoothly  across  the  shock,  as  occurrs  in  most  finite-difference 
solutions,  the  Integration  followed  by  the  differentiation  of  the  term  inside 
the  bracket  can  be  performed  to  give  the  following  approximate  equation. 


(1  -  M*2)(m£  -  M*) 


-  (1 


Y  ~  1 

r  +  1 


m*2)m* 


(60) 
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where  M  =  u/a*  a  $  /a^  and  a*  is  the  sonic  speed  at  M  =  1.  Two 
approximations  have  been  made:  the  (1  -  M*^)  term  on  the  left  side  and  the 

*o  * 

[  1  —  ( y  —  l)/(y  +  1)M  ]M  term  on  the  right  side  are  treated  as  constants 

if 

during  integration,  and  the  relative  change  in  p  is  assumed  to  be  small.  In 
a  fully  isentropic  flow,  p*  is  constant;  the  right  side  of  Equation  (60)  is 
therefore  zero,  and  the  left  side  can  be  rewritten  as  the  familiar  one¬ 
dimensional  potential  equation 

(a2  -  u2)<t»xx  ■  0.  (61) 

A  change  in  p*  occurs  across  a  shock  in  a  real  flow.  The  right  side  of 
Equation  (60)  should  be  related  to  the  error  incurred  at  a  shock  in  a  poten¬ 
tial  flow  calculation.  M*  on  the  right  side  of  Equation  (60)  is  the  average 


it  it  it 

value  of  and  M^,  and  therefore  a  reasonable  approximation  is  M 


1.  If 


this  approximation  is  made,  and  the  left  side  notation  of  Equation  (61)  is  re 
tained  for  clarity. 


la  ‘  '  a* 


/  Ax, 


where  Ax  is  the  axial  spacing  across  the  shock.  Numerical  values  for  the 
right  side  of  Equation  (62)  can  be  obtained  by  introducing  a  Rankine-Hugoniot 
relation. 


Y-l  * 4 

1  -7TTmi 
M*2  -  Hi 

M1  y+1 


it  ft  it 

where  p.  and  p  are  values  of  p  upstream  and  downstream  of  shocks,  respec- 

^  it^ 

tively,  and  is  the  normalized  shock  Mach  number  upstream  of  the  shock. 

This  result  is  used  to  modify  the  conservative  first-  and  second-order 
artificial  viscosities  at  shock  points  as  follows  (see  Equations  (55)  and 


D  , 

J  A 


s  D 


1  -  -T  /  *>,  • 
P 1  / 


where  AX  is  the  distance  between  the  shock  point  and  its  upstream  supersonic 
s 

point  in  the  streamwise  direction. 


4.3  Relaxation  Strategies 

The  finite-difference  approximation  to  Equation  (8)  can  be  solved  by  a 
line- relaxation  scheme  with  the  boundary  conditions  described  in  the  previous 
section.  To  ensure  that  the  relaxation  scheme  corresponds  to  a  convergent 
process,  the  old  and  updated  values  of  the  potential  function,  4  and  4>+,  must 
be  mixed  properly.  The  basic  relaxation  strategies  developed  for  the  present 
method  are  similar  to  the  ones  described  in  References  15  and  18,  except  for 
careful  treatment  of  artificial  viscosities  at  supersonic  and  shock  points. 

In  the  second-order  quasi-conservative  scheme,  the  old  and  new  values  of 
*  contributing  to  the  terms  (q2  -  a2)$sa  +  H  (or  Hs)  of  the  relaxation 
equation  are  chosen  to  ensure  a  convergent  process  in  the  i  =  constant  line 
sweep  according  to 


<q2  -  a2)  »8a  + 


‘  “*.m  ‘  Wi>  ’ ! 


X  (ci,j,k  “  ci-2m,j,k)  +  Rss 
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for  supersonic  points ,  and 


for  shock.  points,  where  c.t,j,k  -  ,  ,  is  the  correction  to  the 

potential  function,  m  is  equal  to  1  or  -  1  if  U  >  0  or  U  <  0,  respectively, 
and  Rgg  is  the  residual  of  the  finite-difference  approximation  to  (q2  -  a2)<f> 
evaluated  using  old  values  of  $.  The  second-order  fully  conservative  scheme 
applies  the  same  strategy  to  treat  the  artificial  viscosities  at  supersonic 
and  shock  points.  A  similar  strategy  also  can  be  developed  for  the  third- 
order  quasi-conservative  and  fully  conservative  schemes. 


5.  ANALYSIS  OF  SHOCKS  IN  CHANNEL  FLOWS 


A  model  problem,  consisting  of  flow  through  a  planar,  symmetric 
converging-diverging  channel  (Figure  2)  was  used  to  evaluate  various  methods 
for  treatment  of  shock  waves  in  potential  flow.  This  configuration  provides  a 
relatively  simple,  inexpensive  framework  for  evaluation  of  numerical  methods, 
prior  to  incorporation  into  a  three-dimensional  method. 

The  flow  becomes  sonic  at  the  throat,  xt,  and  shocks  occur  in  the  diverg¬ 
ing  section  extending  from  xfc  to  xg.  The  slope  of  the  diverging  wall  can  be 
adjusted  so  that  the  predicted  shock  Mach  numbers  range  from  1.0  to  2.0. 

Values  of  pm  vary  from  -0.6  to  0.8,  and  mass  flux  across  each  x  station  is 
computed  to  check  the  mass  flux  conservation. 

A  previously  developed  inlet  program^  was  modified  to  compute  the  flows 
with  shocks  in  the  channel  by  incorporating  the  second-  and  third-order  arti¬ 
ficial  viscosities  and  the  partially  conservative  shock-point  operators  de¬ 
scribed  previously.  A  typical  grid  system  is  presented  in  Figure  3.  In  the 
calculation,  two  meshes  were  used;  the  coarse  mesh  had  50  mesh  elements  in  the 
unwrapped  X-dlrection,  12  mesh  elements  in  the  Burface  normal  direction,  and 
20  equally  spaced  streamwise  mesh  elements  in  the  diverging  section.  The  fine 
mesh  had  twice  as  many  mesh  elements  in  both  directions.  To  ensure  that  a 
converged  solution  in  the  fine  mesh  existed  and  to  improve  the  convergence 


17 


- 1.0  -0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8  2.0 


X 

Figure  3.  Grid  distribution  about  nozzle  C. 
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rate,  the  potential  function  at  the  downstream  outlet  boundary  was  frozen 
after  300  iterations,  while  the  procedure  continued  to  a  total  of  600  itera¬ 
tions.  The  number  of  supersonic  points  usually  ceases  to  change  in  the  last 
50  to  200  iterations.  Solutions  were  obtained  for  various  slopes  of  the 
diverging  wall  and  various  inlet  and  exit  conditions.  After  each  flowfield 
computation,  the  total  mass  flux  across  each  cross-section  was  computed  by 
numerical  Integration: 


where  ^  it,  the  component  of  the  Mach  number  parallel  to  the  nozzle  axis.  The 
nozzle  flow  differs  from  an  external  flow  with  a  shock  in  that  the  flow  down¬ 
stream  of  the  shock  in  the  channel  can  be  regarded  as  a  potential  flow  with 
"to 

p  =  P2  =  constant,  since  the  variation  in  shock  strength  across  the  channel 

is  relatively  small.  For  each  flowfield  computation,  the  mass  flow  rate  and 

a*  were  assumed  constant,  and  the  change  of  the  integral  upstream  and  down- 

* 

stream  of  the  shock  was  interpreted  as  a  change  in  p  across  the  shock. 

Table  1  presents  a  summary  of  results.  Channels  A,  B,  C.  D,  E,  and  F 
have  values  of  diverging  wall  slopes  of  0.15,  0.20,  0.25,  0.40,  0.50,  and 


TABLE  1.  SUMMARY  OF  CHANNEL  FLOW  CALCULATIONS 


0«875,  respectively,  H  represents  the  order  of  artificial  viscosity,  pm  is  the 
parameter,  first  introduced  in  Equation  (55),  controlling  the  degree  of  non¬ 
conservative  differencing,  RH  in  the  pm  column  means  that  the  shock-point 
operators  defined  in  Equations  (64)  or  (65)  were  used,  and  xg  is  the  shock 
location.  Mex^t  is  the  average  Mach  number  at  the  channel  exit;  xe  *  1.495 
for  all  cases  considered  here.  and  M2  are  the  average  Mach  numbers 
upstream  and  downstream  of  the  shock.  By  assuming  M  =>  1  at  the  throat 
(x  -  1),  the  channel  area  ratio  can  be  used  to  find  a  one-dimensional  value  of 
Mj  at  the  shock.  Similarly,  the  area  ratio  between  x  =*  xg  and  x  =  xe  and 
MgXit  can  be  used  to  determine  a  corresponding  value  of  M2.  *-s  M®0*1 

number  downstream  of  the  shock,  given  by  the  Rankine-Hugonlot  relation. 
(p2/p*)c  Is  the  stagnation  density  ratio  across  the  shock,  computed  by  numer¬ 
ically  integrating  the  function  given  in  Equation  (68).  (p^/p*)^  *s  tlie 

analytical  stagnation  density  ratio  across  a  Rankine-Hugoniot  shock  given  by 
Equation  (63),  and  e  is  the  relative  error  of  (p2/p*)c  to  (p^/p^rh*  A®  Pm 
Increases  from  0  to  0.8,  (p2/p*)c  decreases  from  near  unity  to  a  value  less 
than  (p2/p*)rh*  Smaller  values  of  e  generally  mean  better  agreement  between 

M2  and  M2) 

*  RH 

Figure  4  compares  computed  shock  Mach  numbers,  M^  and  M2,  with  Rankine- 
Hugoniot  shocks  and  Mach  numbers  obtained  from  isentropic,  mass-conserving 
relations.  Solutions  obtained  by  setting  Pm  ■  0  always  lie  below  the 
isentropic,  mass-conserving  shock  curve.  Solutions  obtained  by  settisx  pm  * 
0.6  to  0.8  give  reasonable  agreement  with  the  Rankine-Hugoniot  curve  over  a 
wide  range  of  Mj.  Solutions  obtained  by  applying  Equation  (64)  or  (65)  depend 
on  the  determination  of  M*  in  Equation  (63).  However,  in  the  solution 

J. 

process,  M^  is  chosen  to  be  the  largest  value  of  M  at  the  last  two  supersonic 
points  upstream  of  the  shock.  The  value  of  M*  so  chosen  is  always  smaller 
than  the  exact  M*;  because  of  the  smearing  of  the  shock  over  a  few  mesh 
spacings,  therefore,  M^  is  in  general,  overpredicted.  In  case  26,  the  shock 
occurred  at  the  end  of  the  diverging  section,  and  the  prediction  of  M*  agrees 
with  the  Rankine-Hugoniot  value  almost  exactly;  flowfield  gradients  upstream 
of  the  shock  were  small,  resulting  in  a  more  accurate  value  of  M^.  The 
scatter  in  the  computed  solutions  is  believed  to  depend  on  the  variation  of 
Mach  number  across  the  channel,  the  two-dimensionality  of  the  channel  flow, 
and  the  effect  of  mesh  spacing  relative  to  the  shock  orientation. 


Figure  4.  Comparison  of  computed  average  shock 
Mach  number  with  analytical  solutions. 

Figure  5  presents  computed  values  of  p^/p*  compared  with  the  exact 
Rankine-Hugoniot  solution.  The  exact  solution  for  p£/Pi  f°r  an  isentropic, 
mass-conserving  shock  is  unity.  Solutions  obtained  with  pm  *  0  are  all 
slightly  greater  than  unity.  Figures  4  and  5  are  actually  alternate  methods 
of  presenting  the  same  information  because  of  the  unique  relation  between 
and  p^/p*  for  a  particular  value  of  M*.  Figures  6-14  present  tables  and 
line-printer  plots  of  computed  Mach  number  distributions  on  the  wall  and  along 
the  axis  of  symmetry,  and  Rankine-Hugoniot  Mach  number  distributions  for 
various  cases.  The  columns  labeled  MACH-RH  present  the  analytical  one- 
dimensional  Mach  number  distribution;  columns  labeled  MACH-AXIS  and  MACH-WALL 
are  computed  Mach  numbers  distributions  along  the  axis  of  symmetry  and  along 
the  wall,  respectively. 

The  preceding  channel-flow  calculations  demonstrate  that,  in  this 
Instance,  it  is  possible  to  obtain  solutions  from  a  potential  formulation  that 
closely  approximate  solutions  to  the  EuLer  equations.  In  the  flow  about  an 
airfoil  or  wing,  the  shock  strength  is  maximum  near  the  surface  and  decreases 
to  zero  with  increasing  distance  normal  to  the  surface.  Vorticity,  which  is 
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a  R-H  shock-point  operator  , 
Equations  (64)  anJ  (65) 

o  pm  =0.8  partially  conservative 
Apm  =  0  fully  conservative 


— - Rankine-Hugoniot  shocks 

-  lsentropic  mass-conserving  shocks 
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Figure  5.  Comparison  of  computed  stagnation  density 

changes  across  shocks  with  analytical  solutions. 


neglected  in  the  potential  formulation,  is  thereby  introduced  into  the  real 

flow.  Although  the  vorticity  effect  in  transonic  flowfields  is  of  second- 

order,  *  potential  flows  with  variable-strength  shocks  contain  errors  which 

can  be  minimized,  from  an  engineering  standpoint,  but  not  eliminated.  A  good 

on 

example  of  this  is  shown  by  Lock  who  applied  a  partially  conservative  shock- 
point  operator  to  compute  airfoil  flows.  In  the  nonlifting  case,  his 
partially  conservative  solutions  agree  reasonably  well  with  Euler's  solutions, 
while  in  the  lifting  case,  a  discrepancy  persists.  Since  the  present  work 
gives  a  shock-capturing  technique  resulting  in  Mach  number  jumps  that  are 
close  to  the  Rankine-Hugoniot  values,  it  is  possible  to  interpret  the  computed 
velocity  potential  distribution  downstream  of  the  shock  using  the  Rankine- 
Hugoniot  values  of  stagnation  pressure  or  density  at  each  streamline,  rather 
than  with  the  conventional  lsentropic  assumption.  This  interpretation  would 
have  the  effect  of  carrying  the  inherent  error  in  the  solution  from  the 
immediate  vicinity  of  the  shock  wave,  where  mass  and  momentum  are  now  con¬ 
served,  to  the  region  downstream  of  the  shock.  Comparisons  of  potential 
flowfield  computations  with  numerical  solutions  of  the  Euler  equations  are 
needed  to  evaluate  the  usefulness  of  this  interpretation  since  comparisons 
Involving  experimental  data  are  usually  complicated  by  wind-tunnel  wall 


effects  and  viscous-inviscid  interactions. 
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Figure  7.  Mach  number  distribution  of  channel  flow  for  case  7 
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Figure  9.  Mach  number  distribution  of  channel  flow  for  case  13. 
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6.  WING-BODY  FLOWFIELD  COMPUTATIONS 


6.1  Grid  Generation  and  Computational  Domain 

The  finite-difference  approximation  of  the  governing  equation  obtained  in 
the  previous  section  can  be  constructed  with  the  knowledge  of  mesh-point  loca¬ 
tions.  The  coordinate  transformation  derivatives  are  found  at  each  control 
point  within  a  local  second-order  element  shown  in  Figure  1.  Any  scheme  that 
generates  a  grid  system  in  a  regular  computational  domain  can  be  Incorporated 
with  the  finite-difference  equation  solver. 

In  the  present  study,  a  grid-generation  scheme  developed  in  Reference  21 
is  applied.  The  transformation  of  the  physical  space  to  the  computational 
space  is  shown  in  Figure  15.  The  computational  space  is  truncated  at  a  finite 
distance  from  the  wing  surface.  For  the  results  presented  here,  the  far-field 
boundary  is  placed  approximately  five  to  six  root-chord  lengths  from  the  wing 
surface  in  the  streamwise  and  surface  normal  directions,  and  the  spanwise  far 
field  is  located  two  to  three  semi-span  lengths  from  the  wing  tip  and  the 
outboard  far-field.  C-type  meshes  are  generated  which  wrap  around  the  fuse¬ 
lage  nose  and  wing  leading  edge.  Outboard  of  the  wing  tip,  the  mesh  wraps 
around  a  surface  extending  from  the  wing  tip  to  the  outboard  farfield.  De¬ 
tails  of  the  grid-generation  scheme  can  be  found  in  Reference  21. 


(a)  Physical  space  (b)  Computational  space 


GP21-O0O3-16 

Figure  15.  Physical  and  computational  domain  for  a  wing-fuselage  configuration. 
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Typical  grids  used  in  the  present  calculation  are  presented  in  Figures  16 
and  17.  Figure  16  shows  the  grid  distribution  on  an  0NERA-M6  wing  on  a  verti¬ 
cal  wall,  and  Figure  17  shows  the  grid  distribution  on  the  same  wing  on  a 
semi-infinite  fuselage. 

6.2  Boundary  Conditions 

The  necessary  boundary  conditions  include  the  impermeability  condition  on 
the  wing  and  fuselage  surfaces,  the  Kutta  condition  along  the  trailing  edge, 
the  zero  streamwise  variation  on  the  downstream  Trefftz  plane,  and  the  free- 
stream  condition  on  the  other  far field  boundaries.  For  easy  implementation  of 
the  far-field  freestream  condition,  a  reduced  potential,  G,  representing  a 
perturbation  from  the  freestream,  is  introduced  according  to 

♦  (x  cos  a  +  y  sin  a  +  G)  ,  (69) 

where  is  the  freestream  velocity  and  n  is  the  angle  of  attack. 

On  the  boundary  cross-planes,  ACUSA  and  SUNOS,  G  is  set  to  zero,  repre¬ 
senting  the  freestream  condition.  On  the  Trefftz  plane  or  the  boundary  cross¬ 
plane  AGOSA  and  CMNUC,  the  streamwise  variations  are  assumed  to  be  zero; 
therefore,  the  following  two-dimensional  equation  is  applied: 
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Figure  16.  Grid  distribution  on  an  ONERA-M6  Figure  17.  Grid  distribution  on  an  ONERA-M6  wing 
wing  and  a  vertical  wall.  and  a  simplified  fuselage. 


(70) 


C2*YY  +  C3^ZZ  +  C5*YZ  +  Vy  +  Vz 


0  . 


Equation  (70)  is  obtained  from  Equation  (8)  by  neglecting  all  derivatives  in 
the  X-direction. 

On  the  fuselage  and  wing  surfaces,  the  impermeability  condition  is 
applied: 

V  =*  0,  on  the  wing  surface,  (71) 

and 

W  =  0,  on  the  fuselage  surface.  (72) 

Exact  surface-boundary  conditions  can  be  enforced  at  boundary  points  by  sub¬ 
stituting  Equations  (19)  and  (20)  into  Equations  (71)  and  (72),  respectively, 
and  solving  the  equations  for  the  value  of  the  potential  function  at  boundary 
points.  One-sided  differencing  is  used  in  the  surface-normal  direction  so 
that  there  is  no  need  to  extrapolate  the  potential  to  imaginary  points  inside 
the  wing  or  fuselage  surfaces.  However,  Equation  (72)  might  not  be  suitable 
for  highly  distorted  grids  near  the  fuselage  and  wing  intersection.  Boundary 
conditions  obtained  from  the  finite-volume  algorithm1  give  better  results  near 
the  intersection.  Therefore  all  solutions  presented  in  the  subsequent  section 
were  obtained  by  applying  the  finite-volume  surface-boundary  condition  in  the 
cross-plane  ACMG. 

Along  the  trailing  edge,  the  linearized  equation 
(hj  +  h^  +  h*)*^  +  (hj  +  h*  +  hj:)^  +  (h*  +  h*  +  h*)$zz  -  0  (73) 

is  assumed  to  hold.  Equation  (73)  is  obtained  from  Equation  (8)  by  neglecting 
the  nonlinear  velocity  contribution  and  the  cross-  and  first-derivative 
terms.  This  linearized  equation  is  approximately  valid  along  the  trailing 
edge  only  for  wing  cross-sections  having  a  finite  trailing-edge  angle  where 
zero  flow  velocity  can  be  approximately  assumed;  Equation  (73)  can  be  regarded 
as  an  interpolation  operator  when  the  wing  trailing  edge  is  cusped.  The 
circulation  T  at  each  spanwise  location  is  determined  iteratively  as  the 
solution  proceeds.  Constant  discontinuities  in  potential  across  the  cut 
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downstream  of  the  trailing  edge  are  enforced  along  the  sli'eaawise  coordinate 
lines  extending  from  the  trailing  edge  to  the  downstream  far-fleld.  The  value 
of  the  discontinuity  in  each  spanwlse  plane  is  computed  at  the  trailing  edge 
by  satisfying  Equation  (73)  at  both  the  upper  and  lover  trailing  edges. 

Beyond  the  wing  tip,  the  continuity  of  the  potential  function  across  the 
surface  about  which  the  mesh  is  unwrapped  can  be  approximated  by  solving  $Yy  * 
0  at  points  on  this  surface  and  just  next  to  the  tip.  The  same  condition  is 
applied  in  the  FLO-22  code^  to  solve  for  the  potential  function  at  points 
lying  on  the  vortex  sheet. 

6.3  Numerical  Results 

Typical  solutions  obtained  using  the  second-  and  third-order,  quasi¬ 
conservative  and  fully  conservative  schemes  are  presented  in  this  section. 

Two  meshes  are  used  in  all  calculations.  The  coarse  mesh  contains  44  mesh 
cells  in  the  X-direction,  10  mesh  cells  in  the  Y-direction,  and  7  mesh  cells 
in  the  Z-direction,  where  32  x  5  mesh  cells  are  on  the  unwrapped  wing  sur¬ 
face.  The  fine  mesh  has  double  the  number  of  mesh  cells  in  each  direction. 
Two-hundred  relaxation  sweeps  were  performed  on  the  coarse  mesh,  followed  by 
two-hundred  relaxation  sweeps  on  the  fine  mesh. 

Figures  18  and  19  present  comparisons  of  first-  and  second-order  fully 
conservative  solutions  obtained  for  the  0NERA-M6  wing  on  a  semi-infinite 
fuselage  shown  in  Figure  17.  The  fuselage  has  a  constant  radius  of  0.2  semi¬ 
span  length,  measured  from  the  wing  root  to  wing  tip,  in  the  section  of  the 
wing-fuselage  intersection.  There  are  no  experimental  data  for  this  configu¬ 
ration;  however,  experimental  data  are  available  for  the  same  wing  on  a  verti¬ 
cal  wall  in  Reference  23.  Computed  solutions  at  20%  and  65%  semi-span  loca¬ 
tions  are  presented  in  Figures  18  and  19,  respectively,  for  -  0.84  and 
a  -  3.06°,  where  a  is  the  angle  of  attack  and  also  the  angle  of  incidence 
between  the  wing  and  fuselage.  The  fuselage  effect  is  not  pronounced  in  this 
case,  and  agreement  between  the  computed  solutions  and  experiment  is  generally 
good.  At  the  20%  semi-span  location,  the  first-order  solution  agrees  with  the 
second-order  solution  except  for  minor  differences  in  suction  peaks  and  de¬ 
tails  at  the  shocks.  At  the  65%  semi-span  location,  experimental  data  show  a 
distinct  double  shock  on  the  upper  surface.  The  second-order  solution 
obviously  resolves  this  double  shock  better  than  the  first-order  solution, 
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Figure  18.  Comparison  of  first-  and  second-order  fully 
conservative  solutions. 


although  there  are  still  small  discrepancies  between  the  second-order  solution 
and  experiment,  presumably  because  the  mesh  used  for  the  computation  is  rela¬ 
tively  coarse. 

In  Figures  20-24  pressure  distributions  obtained  for  an  0NERA-M6  wing  on 
a  wall.  Figure  16,  are  shown  and  compared  with  experimental  data. The  free- 
stream  Mach  number  is  0.84  and  the  angle  of  attack  is  3.06°.  Second-order 
quasi-conservative  and  first-  and  second-order  fully  conservative  solutions 
were  obtained  at  the  20%,  65%,  and  95%  semi-span  locations  and  are  shown  in 
Figures  20-22.  At  the  20%  semi-span  location,  numerical  solutions  predict 
lower  suction  peaks.  The  plateau  pressures  at  the  20%  and  65%  semi-span 
locations  are  slightly  overpredicted,  while  the  pressures  downstream  of  the 
shocks  are  slightly  higher  than  the  experimental  data.  The  locations  of 
shocks  are  predicted  accurately  by  the  numerical  solutions.  Although  the  mesh 
used  is  still  relatively  coarse,  the  overall  agreement  between  the  numerical 
and  experimental  results  is  satisfactory.  The  quasi-conservative  solutions 


lower  surfaces  of  an  ONERA  wing  on  a 
wall  at  the  20%  semi-span  location. 
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Figure  21.  Pressure  distributions  on  the  upper  and 

lower  surfaces  of  an  ONERA-M6  wing  on 
a  wall  at  the  65%  semi-span  location. 


Figure  22. 


Pressure  distribution  on  the  upper  and 
lower  surfaces  of  an  ONERA-M6  wing  on 
a  wall  at  the  95%  semi-span  location. 
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Figure  23.  Comparison  of  third-order  quasi-conservative 
and  fully  conservative  solutions  at  the  65% 
semi-span  location. 
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Figure  24.  Comparison  of  third-order  quasi-conservative 
and  fully  conservative  solutions  at  the  95% 
semi-span  location. 
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predict  a  more  positive  pressure  at  the  trailing  edge  and  yield  better  agree¬ 
ment  with  experimental  data  upstream  of  the  trailing  edge.  This  result  is 
apparently  attributable  to  enforcement  of  the  exact  surface  boundary  condition 
and  the  linearized  equation,  Equation  (73),  which  approximately  simulates  a 
flow  stagnation  condition  for  finite  trailing-edge  angles,  as  mentioned  pre¬ 
viously.  However,  the  quasi-conservative  solutions  predict  a  more-downstream 
shock  location  at  a  65%  semi-span  location  and  do  not  resolve  the  suction  peak 
as  well  as  the  fully  conservative  scheme.  The  poor  resolution  of  the  suction 
peak  is  caused  by  the  second-order  element  not  resolving  the  surface  curvature 
and  potential  gradient  near  the  leading  edge  as  the  exact  surface-boundary 
condition  Is  enforced.  The  leading-edge  resolution  can  be  significantly 
Improved  by  using  a  third-order  element  as  shown  in  Reference  7. 

Third-order,  quasi-conservative  and  fully  conservative  solutions  are 
obtained  for  the  same  wing  and  compared  with  experiment  in  Figures  23  and  24 
at  the  65%  and  95%  semi-span  locations,  respectively.  The  finite-volume 
boundary  conditions  were  applied  on  the  airfoil  surface  for  both  solutions. 

The  pressure  distributions  in  the  supersonic  region  are  more  accurately 
predicted  at  the  95%  semi-span  location  than  in  the  second-order  solutions 
shown  in  Figure  20,  although  the  plateau  pressure  still  is  underpredicted.  A 
triple  shock  is  found  in  the  computed  pressure  distribution  at  the  65%  semi¬ 
span  location,  although  the  first  shock  is  somewhat  ambiguous.  This  pattern 
also  appears  at  spanwise  locations  between  50%  and  70%.  Evidence  of  the 
triple  shock  is  not  shown  in  the  experimental  data;  however,  triple-shock 
structures  have  been  observed  in  static-pressure  distributions  obtained  with 
other  wings  at  transonic  speeds.  This  feature  probably  is  not  associated  with 
poor  convergence  since  the  solutions  presented  here  have  all  converged  to 
sufficiently  small  residuals.  It  is  also  possible  that  the  triple-shock 
pattern  is  the  result  of  a  numerical  oscillation  caused  by  interaction  between 
the  numerical  frequency  of  the  third-order  artificial  viscosity  and  the  wave 
frequency  associated  with  the  double  shock.  To  explain  the  characteristics  of 
the  third-order  solutions,  finer-mesh  solutions  are  needed. 

A  study  of  the  shock-point  operator  for  the  same  0NERA-M6  wing  on  the 
wall  at  M  -  0.84  and  a  -  3.06  is  shown  in  Figures  25  and  26.  As  mentioned 

flD 

previously,  pm  is  the  parameter  controlling  the  nonconservative  differencing 
of  the  shock-point  operators  in  Equations  (61)  and  (64).  Solutions  are 


obtained  for  pm  -  0,  1,  and  2  at  20%  and  65%  semi-span  locations  by  using  the 
second-order  artificial  viscosity  at  supersonic  points.  A  second-order  fully 
conservative  solution  is  obtained  as  pm  =  0.  As  pm  increases,  the  amount  of 
nonconservative  differencing  increases ,  and  the  additional  mass  flux  intro¬ 
duced  at  the  shock  increases.  By  adjusting  the  value  of  pm,  part  of  shock- 
induced  boundary-layer  displacement  effect  can  be  simulated.  By  increasing 
the  value  of  pm,  the  shock  moves  upstream,  the  shock  strength  decreases,  and 
agreement  of  the  computed  pressure  with  experimental  data  is  significantly 
improved  both  upstream  and  downstream  of  the  shocks.  Solutions  obtained  by 
setting  pm  »  2  seem  to  yield  best  agreement  with  experiments. 
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Figure  25.  Study  of  partially  conservative  shock-point 


operators  at  209o  semi-span  location  on  an 
ONERA-M6  wing. 
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Figure  26.  Study  of  partially  conservative  shock-point 
operators  at  65%  semi-span  location  on  an 
ONERA-M6  wing. 


7.  CONCLUSIONS 


Second-  and  third-order,  fully  conservative  and  quasi-conservative 
schemes  have  been  developed  to  compute  flowfields  about  transonic  wings  and 
wing-body  configurations.  The  quasi-conservative  scheme  was  developed  by 
solving  a  finite-difference  representation  of  a  transformed  full  potential 
equation  formulated  in  this  report  and  enforcing  an  exact  body  surface¬ 
boundary  condition. 

The  second-order  solutions  obtained  have  been  shown  to  provide  better 
resolution  for  a  double  shock  than  the  conventional  first-order  schemes.  The 
third-order  solutions  show  a  triple-shock  pattern.  Additional  study  will  be 
required  to  determine  whether  this  pattern  is  a  real  flowfield  characteristic 
or  a  feature  of  the  numerical  scheme.  The  enforcement  of  an  exact  surface¬ 
boundary  condition  in  the  quasi-conservative  scheme  provides  solutions  with 
better  agreement  with  experiments  upstream  of  the  trailing  edge. 

A  partially  conservative  shock-point  operator  is  introduced  to  control 
the  amount  of  nonconservative  differencing  at  shock  points  and  thus  modify  the 
location  and  strength  of  shocks.  Proper  choice  of  the  shock-point  operator 
significantly  Improves  the  agreement  of  computed  pressure  distribution  with 
experimental  data  near  the  shocks.  A  new  shock  point  operator  is  derived  from 
an  approximate  one-dimensional  flow  analysis  which  properly  considers  the 
stagnation  density  change  across  the  shock  and  predicts  shocks  in  reasonable 
agreement  with  Rankine-Hugoniot  shocks. 
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